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Using local density approximation (LDA) calculations we predict GdFe03-like rotation of Ti06 
octahedra at the n-type interface between LaAlOs and SrTiOs. The narrowing of the Ti d bandwidth 
which results means that for very modest values of U, LDA+?7 calculations predict charge and spin 
ordering at the interface. Recent experimental evidence for magnetic interface ordering may be 
understood in terms of the close proximity of an antiferromagnetic insulating ground state to a 
ferromagnetic metallic excited state. 



PACS numbers: 68.35. Ct, 71.28.+d, 73.20.-r, 75.70. Cn 

Transition metal (TM) oxides in bulk form exhibit 
a huge range of physical properties. Heterostructures 
of TM oxides offer the prospect of greatly enhancing 
these properties or of combining them to realize entirely 
new properties and functionalities. The recent finding 
that a Ti02|LaO interface between the insulating ox- 
ides LaAlOs and SrTi03 can be metallic with an ex- 
tremely high carrier mobility pQ has triggered a surge 
of experimental and theoretical studies of this interface 
[2l[3l[4l[5l[6l[7l[8l[9l HOl HH H2]. The valence mis- 
match at the interface leads to the transfer of half an elec- 
tron per unit cell from LaAlOs (LAO: bandgap 5.6 eV) 
to SrTi0 3 (STO: bandgap 3.2 eV), from the LaO in- 
terface layer to the TiC>2 layer. For a defect-free inter- 
face, the concentration of these "intrinsic" sheet carriers 
is n s heet ~ 3.3 x 10 14 cm -2 . For samples prepared under 
low oxygen pressure (< 10 -4 mbar), the carrier density 
is much larger than this intrinsic value, suggesting their 
properties are determined by extrinsic carriers related to 
oxygen vacancies formed in the STO substrate during the 
growth of the LAO film. 

Samples prepared under increasing oxygen pressure ex- 
hibit a large increase of the sheet resistance [TJ 0J [6] . A 
sheet resistance which decreases on increasing the tem- 
perature from 10K to 70K, large negative magnetore- 
sistance and magnetic hysteresis at low temperatures 
were found in Ref. [4] and these properties are presum- 
ably characteristic of intrinsic interfaces. The partial 
filling of the Ti d states at the n-type interface means 
that locally the electronic structure is intermediate be- 
tween orthorhombic LaTiOs (Ti 3+ ) and cubic SrTiOs 
(Ti 4+ ). LaTi0 3 (LTO) is a G-type antiferromagnetic 
(AFM) Mott insulator with a Neel temperature of 146 K 
P~3j [14j [15] and a band gap of 0.2 eV [16]. Compared to 
cubic STO, the additional electron in the t^g conduction 
band state on the Ti ion leads to a GdFeOs-type crys- 
tal structure (space group Pbnm) which can be derived 
from the ideal perovskite cubic structure by tilting essen- 
tially ideal Ti06 octahedra about the b-axis followed by 
a rotation about the c-axis (see Fig. HI). 



In this paper we suggest that GdFeOs-type distortions 
at LAO | STO interfaces play an essential role in reducing 
the bandwidth of the occupied interface Ti d state, ren- 
dering it very sensitive to on-site Coulomb correlations 
and stabilizing an antiferromagnetic ground state. Sur- 
prisingly little attention has been paid to the possibility 
of such distortions at the LAO | STO interface. Theoret- 
ical studies have been largely concerned with explaining 
the metallic behavior of n-type and the insulating char- 
acter of p-type interfaces [Til 02] H3 HE] HH] with most 
emphasis being placed on the charge transfer between the 
LaO and Ti02 interface layers. 

Though LDA total energy calculations describe equi- 
librium crystal structures with good accuracy and have 
been applied extensively to non-metallic TM oxides [20] , 
we are not aware of applications to the structures of open 
shell correlated solids. To establish a predictive capabil- 
ity for the LAO | STO system, we carried out a system- 
atic study of the crystal structure of the 3d 1 perovskites 
CaVOs, SrVOs, YTi0 3 and LaTi0 3 . We could show 
that LDA calculations reproduce the well-documented 
[2T] GdFeOs-type distortions in this system very well. 
To study various types of magnetic ordering, we took ac- 
count of the Coulomb interaction of the partially filled 
d-shell using a local Coulomb interaction parameter Ud 
[22] which has been shown to describe the ground state 
of YTiOs and LaTiOs qualitatively correctly [23 . The 
LDA+U results are in even closer agreement with exper- 
iment [24] . 

Method: We model the LAO | STO interface using a 
periodically repeated (m, n) supercell containing m layers 
of LAO and n layers of STO. Most of the results to be 
reported below were obtained with the 40 atom c(2 x 2) 
(|, |) supercell depicted in Fig. [I] containing two n-type 
interfaces, and did not change significantly when a (|, |) 
supercell was used instead. To study AFM ordering the 
above supercell was doubled in the xy plane leading to a 
p(2 x 2) unit cell containing 80 atoms. 

The LDA and LDA+Z7 calculations were carried out 
using the projector augmented wave (PAW) method [25] 
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FIG. 1: Bulk LaTi0 3 based LaA10 3 |SrTi0 3 interface struc- 
ture and charge density iso-surface of the surplus electron for 
the charge ordered ferromagnetic state. The orthorhombic 
translation vectors are a, b, c. 



[26] and a plane wave basis as implemented in the Vienna 
Ab initio Simulation Package (VASP) [27j[28]. A kinetic 
energy cutoff of 500 eV was used and the Brillouin zone 
of the 40 atom supercell was sampled with a 6x6x4 k- 
point grid in combination with the tetrahedron method 
[29] . Compared to the LDA, the LDA+/7 approach gives 
an improved description of d electron localization [22] . 
We use the rotationally invariant LDA+Z7 method [30] 
with, unless stated otherwise, Ud = 5 and Jd = 0.64 eV 
for Ti d states. The La / states are forced to lie higher in 
energy by imposing a large Uf — 11 eV and Jf = 0.68 eV 

ED- 

Results: Because samples are grown on STO sub- 
strates, we fix the in-plane lattice constant at the ex- 
perimental value for STO. This imposes a strain on LAO 
with a 3% smaller bulk lattice constant. Our starting 
point is a fully relaxed, tetragonal, non-magnetic metallic 
(RT-NMM) p(l x 1) configuration. We then remove sym- 
metry constraints to obtain, in order of increasing energy 
gam (Fig. |2Jl): a c(2 x 2) distorted non-magnetic metal- 
lic (D-NMM) structure; a metastable spin-polarized, non 
charge-ordered (NCO) state, with a single interface-Ti 
charge state, obtained by switching on Uj z ; a charge- 
ordered ferromagnetic insulating (CO-FMI) state with an 
energy gap of 0.44 eV in which the Ti sites are not equiva- 
lent (Fig.[T]); ap(2x2) charge-ordered anti- ferromagnetic 
insulating state (CO-AFMI). 

The charge transferred from LAO to the smaller 
bandgap STO is bound close to the interface. As pre- 
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FIG. 2: (a) Energy gained per surplus electron by relaxing 
symmetry constraints. The LDA+?7 energies (dashed lines) 
and energy differences calculated here with Ud = 5 eV depend 
on the value of U used and so are just indicative, (b) Phase di- 
agram as function of value of Uj l . Top: energy gap. Middle: 
disproportionation order parameter rj = \m — 7i2|/(ni + 712) 
where rti is the occupation of the Ti d states on ion i. Bottom: 
AE = \Ecofm—Ecoafm\- (c) Spatial magnetic distribution 
of the COAFM states {Ud = 4eV) at the interface. 



dieted by Okamoto and Millis [17] and confirmed for 
LTO|STO interfaces using calculations similar to those 
presented here [31], the effect of lattice relaxation is to 
smooth the potential discontinuity and allow charge to 
leak into the bulklike STO. For our NCO-FMM state, 
the occupied part of the bottom of the conduction band 
is almost leV wide and we confirm the occurrence of 
charge leakage at the LAO | STO interface. 

The NCO-FMM state is unstable with respect to 
charge disproportionation into inequivalent Ti 3,2+ and 
Ti 38+ interface ions as a result of the Coulomb repul- 
sion and spin polarization on the Ti ions favoring integer 
occupancy; the transferred electron is localized on the 
nominally 3+ interface Ti 3+ ions with little leakage into 
bulk STO (Fig. [3|. The atoms relax to accommodate 
the charge ordering: small displacements of the oxygen 
atoms surrounding the Ti 3+ ion lead to that Ti06 octa- 
hedron expanding while the other one contracts. The 
combination of this distortion with the potential bar- 
rier and reduced symmetry at the interface leads to one 
of the Ti t2 g states splitting off completely from the 
bottom of the STO conduction band and a gap open- 
ing up for quarter filling. The transferred electron re- 
mains strongly localized at the interface with essentially 
no leakage into the bulk layer. Because of the GdFeOs- 
type interface distortion, the occupied Ti state has 
Q.79\xy) + 0.33\yz) +0.52|xz) orbital character (Fig.0. 

In a p{2 x 2) geometry, the kinetic energy can be fur- 
ther reduced by flipping the spins of equivalent Ti ions 
in a checkerboard pattern. The bandwidth is narrowed 
to 0.4 eV and the energy gap of the CO-AFMI state in- 
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FIG. 3: Interface densities of states (DoS) for the lowest en- 
ergy charge ordered antiferromagnetic insulating (CO-AFMI) 
configuration, (a) and (b) projected 3d DoS of (nominally) 
Ti 3+ and Ti 4+ interface Ti ions. The transferred electron has 
tig (red) character with the e g states (black) far above the 
Fermi energy; (c) projected 3d DoS of bulk Ti 4+ ions; (d) 
total DoS. 

creased by 20 meV compared to the CO-FMI state. The 
total energy is lowered but the energy gain ~ 10 meV is 
small, close to the limit of our accuracy. 

The effect of varying Uj z between 2 and 7 eV, keep- 
ing J fixed at 0.64 eV, is shown in Fig. [2}d. Between 
3.5 and 4 eV, charge disproportionation takes place, AF 
ordering becomes more favourable and a metal insulator 
transition occurs. Increasing Ud leads to a larger CO gap 
and with increasing charge disproportionation the stripe 
charge and AF ordered pattern shown in Fig. evolves 
into a checkerboard pattern. The energy difference be- 
tween CO-FMI and CO-AFMI states reaches a maximum 
at Ud = 5eV, the effect of increasing U being to suppress 
the exchange coupling between neighbouring spins. The 
most important result of this study is that the Ti06 rota- 
tion reduces the bandwidth of the split-off band so that 
very modest values of Ud result in FM and AFM order- 
ing and a metal insulator (MI) transition close to the 
FM-AFM phase boundary. 

Structure: Because the LAO lattice constant is 3% 
smaller than that of STO, LAO might be expected to 
shrink in the z direction when it is forced to match in- 
plane to STO. However, at the interface, we find that 
the La and Sr planes are separated by 4. 003 A which is 
substantially larger than (3.905 + 3.789)/2, the arith- 
metic average of the LAO and STO lattice constants, 
and somewhat larger than the lattice constant of a fic- 
titious cubic LTO phase with the same volume as the 
experimentally observed orthorhombic phase, but con- 
sistent with recent experiments [32]. By comparing re- 
sults for (§, |) and (§, \) multilayers for the NCO-FMM 
and CO-FMI states, we find that the GdFeOs rotation 
of the TiOe octahedra is confined to the interface. The 
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FIG. 4: Sketch of the vertical displacement patterns of the 
Ti and O interface atoms for (a) an ideal cubic structure (b) 
a relaxed p(l x 1) structure (c) an ideal c(2 x 2) "GdFeOs" 
structure and (d) the fully relaxed c(2 x 2) "buckling" struc- 
ture. 

central STO (LAO) layers in the (|, |) multilayers are 
essentially cubic (tetragonal). Compared to a bulk or- 
thorhombic LaTiOs with an appropriately scaled lattice 
parameter, the shift of the La and Sr ions is reduced, 
presumably as a result of the constraints imposed by | 
or I layers of cubic STO and the availability of only half 
an extra electron in the Ti d states. 

Previous LDA and LDA+Z7 calculations found 
ferroelectric-like displacements of negatively charged O 
and positively charged Ti ions when the geometry was 
optimized for p(l x 1) structures [I2j [31] resulting in a 
zigzag pattern of O-Ti-O-Ti-0 Ti displacements perpen- 
dicular to the interface; see Fig. [4)3. Combining this with 
the GdFeOs rotation of the TiOe octahedra (Fig. [4^) 
leads to the buckling pattern shown in Fig. |4]i in which 
one oxygen is displaced by ~ 0.15 A out of an otherwise 
essentially linear chain. 

Discussion: Let us assume that we are just on the insu- 
lating side of the MI transition in Fig. i.e., Uj l ~ 4 eV. 
Our calculations predict that the strong correlation of Ti 
d states will lead to electrons transferred from the LaO 
to the Ti02 layer being trapped on the interface Ti ions 
in an AFMI state. The localization of these intrinsic car- 
riers would explain the high sheet resistance (> 10 4 ^/D) 
and low carrier concentration (~ 10 13 cm -2 ) of samples 
prepared at high oxygen pressure [TJ |H |6j [7] . The exis- 
tence of a band gap would explain the decrease of sheet 
resistance with increasing temperature [4 . Application 
of an external magnetic field would force the AFMI sys- 
tem into a ferromagnetic insulating state with a smaller 
band gap, leading to a reduced resistance as observed. 

This cannot be the whole picture, however. For exam- 
ple, a simple (doped) semiconductor model would predict 
carrier freeze-out at low temperatures and this is not ob- 
served. The role of oxygen vacancies as well as that of 
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the cation mixing at the interface needs to be clarified. 
If we focus on the ideal of an intrinsic, abrupt interface, 
perhaps the most pressing issue to be resolved experi- 
mentally relates to the appropriate value of Uj l which is 
a free parameter in our otherwise parameter-free study. 
A spectroscopic study identifying the existence and size 
of the predicted CO gap would be invaluable for making 
further progress. 

Summary and Conclusions: Using LDA and LDA+/7 
calculations, we have shown that the half electron (per 
interface Ti ion) transferred from the LaO layer on one 
side of an LAO|STO interface to the Ti02 layer on the 
other side favors a rotation of Ti06 octahedra just as it 
does in bulk LTO. As in bulk LTO, the distortion is cru- 
cial to the formation of a charge-ordered AFMI ground 
state and other charge, magnetic and orbital properties, 
and results in a characteristic buckling of the Ti-O-Ti 
bonding at the interface. Ionic relaxation plays a crucial 
role in determining the localization of the extra electron 
on the interface Ti ions and leakage into the bulk layer. 
However, charge ordering suppresses this leakage even 
when relaxation is included and, in the CO states, the 
electrons are strongly localized at the interface. Our cal- 
culations suggest an explanation for recent experimental 
results in terms of the proximity of an AFMI ground 
state to a FMM excited state (or a FMI state with re- 
duced band gap). With only half an electron trapped in 
an otherwise empty Ti d orbital, the LAO|STO interface 
is an attractive object for quarter-filled Hubbard model 
studies. 
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